Field-induced long-range magnetic order in the spin-singlet ground state system 

YbAlsCs: Neutron diffraction study 
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The 4/-electron system YbAlsCs with a non-magnetic spin-dimer ground state has been studied 
by neutron diffraction in an apphed magnetic field. A long-range magnetic order involving both 
ferromagnetic and antiferromagnetic components has been revealed above the critical field He ^ 6T 
at temperature T=0.05K. The magnetic structure indicates that the geometrical frustration of the 
prototype hexagonal lattice is not fully relieved in the low-temperature orthorhombic phase. The 
suppression of magnetic ordering by the remanent frustration is the key factor stabilizing the non- 
magnetic singlet ground state in zero field. Temperature dependent measurements in the applied 
field H=12T revealed that the long-range ordering persists up to temperatures significantly higher 
than the spin gap indicating that this phase is not directly related to the singlet-triplet excitation. 
Combining our neutron diffraction results with the previously published phase diagram, we support 
the existence of an intermediate disordered phase as the first excitation from the non-magnetic 
singlet ground state. Based on our results, we propose YbAlsCs as a new material for studying the 
quantum phase transitions of heavy- fermion metals under the influence of geometrical frustration. 

PACS numbers: 75.25.-j 



The formation of a non-magnetic singlet ground state 
due to spin dimerization is an extremely rare phe- 
nomenon in 4/-electron systems. It requires strong quan- 
tum effects which are usually significant only in low di- 
mensional spin = 1/2 systems. Due to a larger total 
angular momentum J, and the three-dimensional char- 
acter of interactions via conduction electrons, the spin 
dimer state is not favourable in most rare earth com- 
pounds. The exceptional cases are Yb4As3 where an 
effective spin S = 1/2 Heisenberg chain is believed to 
be realized [Ij] and the recently proposed spin-dimer sys- 
tem YbAlsCs [l-lsl- Due to the effect of the crystalline 
electric field, the ground state Kramers doublet is sug- 
gested to be well isolated from the exited states in these 
compounds and the 4/ electrons at low temperature are 
expected to behave as a 5 = 1/2 spin system. 

The non-magnetic nature of the ground state in 
YbAlsCs was initially interpreted as an antiferro- 
quadrupolar ordered state [6] which takes place at 
T^ ~80K, where the specific heat exhibits a A- type 
anomaly. Later, this interpretation was discarded by 
Ochiai et al. [2] who observed a similar sharp peak in the 
specific heat for the Lu-based analogue at 110 K and at- 
tributed it to a structural phase transition present in both 
compounds. The new interpretation proposed by Ochiai 
et al. [2] implies formation of isolated dimers due to the 
structural distortions which promote the non-magnetic 
spin-singlet ground state. This idea explains the low tem- 
perature specific heat and magnetization data assuming 
the spin gap to be ~15K. More direct evidence of the 



singlet-triplet excitations in YbAlsCs has been presented 
by Kato et al. [3|] and Adroja et al. [4] based on inelastic 
neutron scattering. 

By analogy with insulating d-electron dimer systems, 
one can expect field-induced long-range magnetic order in 
YbAlsCs at the critical magnetic field closing the spin- 
gap due to the Zeeman effect. It has been shown that 
this kind of quantum phase transition can be modeled 
as a Bose-Einstein condensation in a system of weakly 
interacting bosons 0, HI. At the same time, the metal- 
lic nature of YbAlsCs introduces an interplay between 
RKKY and Kondo interactions which, in turn, suggests 
that the magnetic field will induce behavior that is dis- 
tinct from that in insulating dimer systems. The aim of 
the present study, therefore, was to directly seek and in- 
vestigate the field-induced long-range magnetic order in 
YbAlsCs through neutron diffraction. Our data indeed 
reveal a change of the ground state from non-magnetic 
below the critical field Hc=6T to magnetically ordered 
above He. Surprisingly, the suggested magnetic struc- 
ture does not imply a strong exchange coupling for the 
expected isolated dimers and indicates the presence of 
geometrical frustration in spite of the lack of the three- 
fold symmetry in the low temperature structural phase. 
Based on the obtained results, one can conclude that the 
remanent frustration is the crucial ingredient promoting 
the low-temperature singlet ground state and is respon- 
sible for the unusual excitations in the system; in turn, 
our results suggest that YbAlsCs provides a new setting 
to study heavy-fermion quantum phase transitions under 
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FIG. 1. (Color online) (a) Orthorhombic structure of 
YbAlsCs (the parent hexagonal cell is shown as well). The 
primary orthorhombic displacements of Al and C are repre- 
sented by arrows, (b) Neutron diffraction patterns collected 
above and below the structural phase transition, (c) Ri- 
etveld refinement of the neutron diffraction pattern collected 
at T=2K. Inset shows the unit cell parameters as a func- 
tion of temperature, (d) Triangular layer formed by Yb^^ 
ions and the exchange parameters in the orthorhombic phase. 
The interatomic Yb-Yb distances (evaluated at T=2K) cor- 
responding to the non-equivalent exchange parameters are: 
3.3970(1)A for Ji{i = 1-3) and 3.3995(1) for J4. J2-bonds 
are arbitrarily taken to show the isolated dimers as ellipses 
around the bonds. 



the influence of geometrical frustration. 

The high resolution neutron diffraction experiments 
were performed on a 2g polycrystalline sample YbAlaCa 
prepared as described by Kosaka et al. [6] The mea- 
surements were carried out on HRPD (in zero mag- 
netic field and 2K<T<300K temperature range) and 
the WISH diffractometer |9] (in the 0T<H<12T field 
range at T=0.05K and in the 0.05T<T<20K tempera- 
ture range at H=12T) at the ISIS Facility of the Ruther- 
ford Appleton Laboratory (UK). To avoid alignment of 
sample grains, in the field-dependent measurements, the 
powder was cold pressed into pellets fitting tightly inside 
a 6mm vanadium can. The refinement of the crystal and 
magnetic structures were done using the FullProf soft- 
ware [iq| . 

The high temperature neutron diffraction patterns 
(T>80K) can be successfully refined in the hexago- 
nal P6s/mmc space group consistent with the previous 
diffraction studies p, luj . The structural parameters ob- 
tained at T=100K are summarised in Table III The Yb^+ 
ions create two-dimensional triangular layers (Fig. [TJi) 
causing the geometrical frustration for the associated 4/ 
magnetic moments. The low dimensionality of the sys- 



TABLE I. Structural parameters of YbAlsCs refined in 
the hexagonal {PGs/mmc, a=3.39804(l)A, c=17.12430(9)A, 
RBragg=4..18 %) and orthorhombic {Pbca, a=17.11895(9)A, 
6=6.79900(10)A, c=5.88243(8)A, RBragg=4.A9 %) phases. 
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All 
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1.65(5) 


CI 


4/ 


1/3 
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1.19(2) 


C2 


2c 


1/3 


2/3 
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1.19(2) 
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3/4 


0.85(2) 


Alll 
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0.13297(43) 


1/8 
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0.92(3) 


A112 
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0.13159(43) 
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0.92(3) 


A12 
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1/8 
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8a 
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1.03(2) 
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8a 


0.59529(13) 
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11/12 


1.03(2) 


C2 


8a 


0.25023(28) 


1/8 


11/12 


0.96(2) 



tem is ensured by the large inter-layer distance ~ 8. 6 A 
resulting in a significant predominance of the in-plane in- 
teractions over the out-of-plane ones. Below T^, a set of 
very weak additional reflections appears (Fig. [Ud), indi- 
cating the structural phase transition. In agreement with 
the single crystal X-ray diffraction work by Matsumura et 
al. [12], these reflections can be indexed with one of the 
symmetry related propagation vectors: ki = 1/4, 1/4,0, 
k2 = -1/4,1/2,0 or ks = -1/2,1/4,0. However, the 
atomic coordinates reported by these authors did not 
work with our diffraction data. To get a better structural 
model a detailed symmetry analysis has been performed. 
The above mentioned wave vector star combines four ir- 
reducible representations, Ai(z = 1 — 4). For each repre- 
sentation, isotropy subgroups and displacive modes were 
generated using ISOTROPY [13] and ISODISTORT [Til 
software and checked in the refinement procedure versus 
the experimentally measured superstructure reflections. 
We restricted our analysis by considering only subgroups 
associated with a single k and — k pair as experimentally 
found by Matsumura et al. jl2[ . The analysis resulted in 
the conclusion that the displacive modes involving only 
Al and C (mainly A12 and CI) and associated with the 
A3 representation and the (a, 0, 0, a, 0, 0) order parameter 
direction drive the phase transition at T^. These primary 
modes have the Phca symmetry in agreement with the re- 
flection conditions deduced by Matsumura et al. [12] and 
induce the atomic displacements along the hexagonal c- 
axis (Fig. [1^). We found the Yb ions were not displaced 
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FIG. 2. (Color online) (a) Neutron diffraction patterns col- 
lected at different temperatures and magnetic fields. Arrows 
indicate positions of the additional reflections induced by the 
magnetic field, (b) Integrated intensity of the antiferromag- 
netic peak (210) as a function of the magnetic field, (c) Fer- 
romagnetic contribution to the (200) nuclear reflection as a 
function of the magnetic field, (d) Refinement of the magnetic 
intensities obtained by subtraction of the patterns collected 
at T=0.05K in H=12T and H=OT (magnetic Bragg factor 

-Rsra^^ = 6.2%) . 



within the precision of our diffraction experiment. Tak- 
ing into account this result the final refinement (Fig. [ih) 
was done in the orthorhombic Pbca space group (related 
to the hexagonal PQ^/mmc by: ao = Ch, bo = 2ah + 2bh 
and Co = bh — ah) using only four parameters varying 
the X coordinates for the six nominally independent Al 
and C sites. The y and z coordinates for these atoms 
as well as all coordinates for Yb were fixed to their high 
symmetry values corresponding to the high-temperature 
hexagonal phase (Table H]) . This simplified approach to 
the refinement of the low-temperature phase is the only 
one possible due to the limited number of superstructure 
reflections in the powder data allowing us to determine 
only primary distortions. 

Considering the effect of the symmetry lowering 
on the exchange topology of the Yb sublattice, the 
favourable symmetry conditions to form isolated dimers 
should be pointed out. In the low temperature or- 
thorhombic phase there are four non-equivalent in-plane 
exchange couplings shown in Figure [T]i. J 2 and J3 form 
isolated dimers and can potentially promote the non- 
magnetic singlet ground state. The corresponding Yb- 
Yb inter-atomic distances vary only slightly due to the 
small variations of the a and h unit cell parameters (Fig. 
[It inset), which indicates that the main factor renormal- 
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FIG. 3. (Color online) (Top) Two degenerate spin configu- 
rations denoted as A and B in a triangular layer. (Bottom) 
Superposition of the degenerate configurations taking with 
same (left) and opposite signs (right). 



izing the exchange parameters should be related to the 
Al and C displacements. 

The absence of any sign of long-range magnetic order- 
ing in our diffraction data down to 0.05K is fully consis- 
tent with the previous neutron diffraction , Mossbauer 
spectroscopy [15| and muon spin relaxation [16] studies. 
However, application of magnetic field ~ 6T changes the 
non-magnetic ground state of YbAlaCa. A clear ferro- 
magnetic signal on top of some nuclear peaks and a set 
of additional reflections in a low-Q region (Fig. [2^) in- 
dicate the onset of long-range magnetic order. The field 
dependence of both ferromagnetic (Fig. [2t) and antifer- 
romagnetic (Fig. ^) components demonstrates a crit- 
ical behaviour with the critical exponents 0.24(3) and 
0.32(6), respectively. The quantitative refinement of the 
magnetic intensity (Fig. [2]i) revealed a uniform ferro- 
magnetic component along the orthorhombic c-axis and 
an antiferromagnetic component along the 6-axis. The 
latter does not require enlargement of the orthorhombic 
nuclear unit cell and all magnetic peaks can be indexed 
with the k = propagation vector. Using representation 
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FIG. 4. (Color online) Two examples of the field- induced 
magnetic structures combining degenerate A and B triangu- 
lar layers. The monoclinic P2[ structure (a) is fully ordered, 
the ferromagnetic (along the c-axis) and antiferromagnetic 
(along the 6- axis) components equal 1.6(2) fiB and 1.6(2)/is, 
respectively, in the magnetic field H=12T. The orthorhombic 
P2i2i2i structure (b) consists of fully ordered ferromagnetic 
component, 1.6(2)/iB, and partially disordered antiferromag- 
netic one, 2.2 (2) /is (only half of the sites carry non-zero an- 
tiferromagnetic component). 



theory we classified different magnetic configura- 

tions according to irreducible representations of the Pbca 
space group and used them in the refinement procedure. 
No solution within a single irreducible representation has 
been found, therefore combinations of different represen- 
tations were tested. We found that many combinations 
describe the observed magnetic intensities equally well 
(Fig. [2li) and no unique solution can be deduced from 
the powder data. However, the important point is that 
all these solutions can be presented as admixtures of two 
configurations denoted as A and B in Figure [3] (top) 
and representing magnetic ordering in the triangular lay- 
ers. It can be a simple interchanging of these layers, 
[AJa^^o [^]a:~o.5 ^ ^ ' • ' aloug the orthorhom- 

bic a-axis (the former hexagonal c-axis), which in combi- 
nation with the ferromagnetic component results in the 
canted structure (Fig. |4^) with monoclinic P2[ symme- 
try. Or a more complex orthorhombic P2'^2'^2i combina- 
tion involving disordered antiferromagnetic components 
on half of the Yb sites (Fig. [3] bottom) and implying 
[A - B]a:r^o ^ [A + B]^^o.5 ^ " ^ ' ' ' alterna- 

tion (Fig. HId). Practically, it means interchanging two 
types of the triangular layers; one of them is randomly 
represented by the A or B configurations taken with dif- 
ferent signs, the second one is randomly represented by 



these configurations taken with same sign. Note that 
the [B - A]^^o ^ [A + B]^^o.5 ^ [B - A]^^i ^ • • • 
stacking type is not equivalent to the previous one and 
has a distinct magnetic symmetry but reproduces equally 
well the diffraction pattern. There are many intermediate 
situations with sites half fully-ordered and half partially- 
ordered (intermediate between the structures shown in 
Figured^ and[4)3). However, the most symmetric variant 
(Fig. [4)3) is considered to be preferable since it involves 
only two representations and therefore a lower degree of 
degeneracy. The value of the ordered moments for Yb^+ 
ions depends on the model and varies from 2.3(3) jj^b for 
the fully-ordered configurations up to 2.7(3) fiB for the 
partially-disordered ones with only half of the sites car- 
rying the ordered antiferromagnetic component. 

The admixture of the A and B configurations signifies 
their degeneracy and points to the fact that the geomet- 
rical frustration is not fully relieved by the structural 
phase transition at T^. Quantitatively, the frustration 
can be expressed by the ratio ( J2 + -^3)/2 = Ji between 
the in-plane exchange parameters, which is held in the or- 
thorhombic phase. This ratio is not consistent with the 
strong predominance of either J2 or J3 exchange coupling 
expected from the symmetry arguments to form isolated 
dimers. A possible explanation of this experimental fact 
is that the ground state of Yb^+ ions in the field- induced 
ordered phase is essentially different from the zero-field 
singlet state. This assumption comes from the large or- 
dered moment of Yb obtained from the experimentally 
measured magnetic intensities. The new ground state 
adopts the bigger moment to gain full advantage of the 
magnetically ordered state and can renormalize the ex- 
change parameters in the system. 

The results obtained indicate directly that the frus- 
trated nature of the orthorhombic phase is the key fac- 
tor to stabilize the non-magnetic singlet ground state in 
YbAlsCa as it has been originally suggested by Kato et 
al. [3] The magnetic fluctuations between degenerate 
manifolds caused by the frustration prevents the system 
from choosing a unique ordered pattern and the spin- 
dimerization takes place to lift the degeneracy of the 
ground state. 

Field-induced transitions to magnetically ordered 
phases are commonly observed in many d-electron dimer 
systems [8]. However, in strong contrast with these sys- 
tems, YbAlaCa does not show a cusp in magnetization 
curves measured in fields H>Hc as a function of tem- 
perature [sf, which can be attributed to the onset of 
magnetic ordering at the critical temperature TAr(H). 
To clarify this important issue, we performed tempera- 
ture dependent diffraction measurements in the magnetic 
field H=12T applied at 0.05K. Surprisingly, both the fer- 
romagnetic and antiferromagnetic contributions to the 
Bragg peaks were observed up to the highest measured 
temperature T=20K (Fig. [2^). A detailed temperature 
study in different magnetic fields is in progress and will be 
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reported elsewhere. This observation indicates that the 
field-induced long-range magnetic order is not directly re- 
lated to the singlet-triplet excitation and can be induced 
at temperatures significantly higher than the spin gap in 
agreement with the magnetization data recently reported 
by Hara et al. [5] 

Taking into account the frustrated nature of the 
YbAlaCs crystal structure and the strong antiferromag- 
netic interactions (the paramagnetic Curie temperature 
is ~-100K), the system is expected to be in a "classical" 
spin-liquid regime in a wide temperature range below T^. 
Thus, the field-induced long-range magnetic order can 
be associated with the properties of this phase rather 
than with the singlet-triplet excitation. This considera- 
tion presumes that the spin dimerization takes place from 
the disordered spin-liquid phase, in other words, the non- 
magnetic quantum singlet state competes with this phase 
and therefore, one can expect that the first excitation, 
when the spin gap is closed, is the disordered spin-liquid 
phase as well. The long-range ordering therefore is the 
second excited state. The presence of the intermediate 
field-induced disordered phase at low temperatures can- 
not be directly deduced from the powder diffraction data 
since this phase does not contribute to the Bragg diffrac- 
tion, but the recent single crystal magnetization data re- 
ported by Hara et al. [5] clearly indicates the presence 
of the intermediate phase below 6T. This intermediate 
phase indeed has been recently confirmed through inelas- 
tic neutron scattering study in applied magnetic field by 
Adroja et al. Q, which revealed that the position and in- 
tensity of the inelastic excitations changed dramatically 
between H=3 and 5T. 

The metallic nature of YbAlaCs suggests the impor- 
tance of Kondo effect and its interplay with RKKY in- 
teractions. We therefore expect it to be a new material 
that can be used to study the global phase diagram of 
antiferromagnetic heavy- fermion metals [{3, [HI- There 
is a burst of recent interest in studying such a phase 
diagram through materials that can tune the degree of 
local-moment quantum fiuctuations through dimension- 
ality pJi] or geometrical frustration [20-22] . YbAlsCa be- 
longs to the latter class of materials, with the distinction 
that the relevance of geometrical frustration is explicitly 
established even though the system is metallic. Further 
support for our picture comes from the observation of a 
large value of the Sommerfeld coefficient (7 = C /T) [5] 
at intermediate magnetic fields. 

In conclusion, the geometrically frustrated system 
YbAlaCa exhibits a first order structural phase transition 
from hexagonal P63 / mmc to orthorhombic Phca symme- 
try at T^ ~80K. The primary order parameter driving 
the transition involves displacements of Al and C ions 
along the hexagonal c-axis. The structural distortions 
renormalize the exchange parameters promoting the for- 
mation of isolated dimers but the orthorhombic symme- 
try does not fully release the frustration. Suppression of 



magnetic ordering by the remanent frustration stabilizes 
the non-magnetic singlet ground state. Application of 
a magnetic field above He ^6T induces long-range mag- 
netic order at T=0.05K. The magnetic structure involves 
a homogeneous ferromagnetic component along the or- 
thorhombic c-axis and an antiferromagnetic component 
along the 6-axis. The latter is likely to be disordered 
on half of the Yb sites. In the magnetic field H=12T, 
both the ferromagnetic and antiferromagnetic compo- 
nents persist up to 20K indicating that the long-range or- 
der is not directly related to the singlet-triplet excitation. 
A field-induced intermediate disordered phase is likely to 
exist as the first excitation from the non-magnetic singlet 
ground state. This opens the primary question whether 
the long-range magnetically ordered phase in YbAlaCa 
really displays similar physics to that of a Bose-Einstein 
condensation of magnons and how the Kondo effect in- 
terplays with the magnetic frustration. 

Acknowledgment: We would like to thank A.D. Hillier, 
J.-G. Park, J.R. Stewart, T. Guidi and P. Riseborough 
for interesting discussions. DTA would like to acknowl- 
edge financial assistance from CMPC-STFC Grant No. 
CMPC-09108. The work at Rice has been supported 
in part by the NSF Grant No. DMR-1006985 and the 
Robert A. Welch Foundation Grant No. C-1411. 



* email: dmit ry. khalyavin@st f c . ac . uk| 
^ email : de vashibhai . adr o j a@ st f c . ac . uk] 
[1] M. Kohgi, K. Iwasa, J.-M. Mignot, A. Ochiai, and T. 

Suzuki, Phys. Rev. B 56, R11388 (1997). 
[2] A. Ochiai, T. Inukai, T. Matsumura, A. Oyamada, and 

K. Katoh, J. Phys. Soc. Jpn. 76, 123703 (2007). 
[3] Y. Kato, M. Kosaka, H. Nowatari, Y. Saiga, A. Yamada, 

T. Kobiyama, S. Katano, K. Ohyama, H. S. Suzuki, N. 

Aso, and K. Iwasa, J. Phys. Soc. Jpn. 77, 053701 (2008). 
[4] D.T. Adroja et al, ISIS Experimental Reports RB920466 

(2010) and RB1210320 (2013). 
[5] K. Hara, S. Matsuda, E. Matsuoka, K. Tanigaki, A. 

Ochiai, S. Nakamura T. Nojima, and K. Katoh, Phys. 

Rev. B 85, 144416 (2012) 
[6] M. Kosaka, Y. Kato, C. Araki, N. Mori, Y. Nakanishi, 

M. Yoshizawa, K. Ohoyama, C. Martin, and S. W. K. 

Tozer, J. Phys. Soc. Jpn. 74, 2413 (2005). 
[7] T. Nikuni, M. Oshikawa, A. Oosawa, and H. Tanaka, 

Phys. Rev. Lett. 84, 5868 (2000). 
[8] T. Giamarchi, C. Rueegg, and O. Tchernyshyov, Nature 

Physics 4, 198 (2008). 
[9] L.C. Chapon, P. Manuel, P.G. Radaelli, et al., Neutron 

News 22 22 (2011). 
[10] J. Rodriguez Carvajal, Physica B 193, 55 (1993). 
[11] T. M. Gesing, R. Potgen, W. Jeitschko and U. Wort- 

mann, J. Alloys Compd. 186, 321 (1992). 
[12] T. Matsumura, T. Inami, M. Kosaka, Y. Kato, T. Inukai, 

A. Ochiai, H. Nakao, Y. Murakami, S. Katano, and H. 

S. Suzuki, J. Phys. Soc. Jpn. 77, 103601 (2008). 
[13] H. T. Stokes, D. M. Hatch, and B. J. Campbell, isotropy, 

[stokes.byu.edu/isotropy.html] . 



6 



[14] B. J. Campbell, H. T. Stokes, D. E. Tanner, and D. M. 

Hatch, J. Appl. Crystallogr. 39, 607 (2006). 
[15] S. K. Dhar, P. Manfrinetti, M. L. Fornasini, and P. 

Bonville, Eur. Phys. J. B 63, 187 (2008). 
[16] T. Yoshida, T. Endo, H. Fujita, H. Nowatari, Y. Kato, M. 

Kosaka, K. Satoh, T. U. Ito, and W. Higemoto: Abstr. 

Meet. Physical Society of Japan (2006 Autumn Meet.), 

Part 3, p. 414, 23pPSA-34. 
[17] Q. Si, Physica B 378, 23-27 (2006); Q. Si, Phys. Status 

Solidi B 247, 476-484 (2010). 



[18] P. Coleman and A. H. Nevidomskyy, J. Low Temp. Phys. 

161, 182-202 (2010). 
[19] J. Custers et al. Nature Mater. 11, 189 (2012). 
[20] M. S. Kim and M. C. Aronson, Phys. Rev. Lett. 110, 

017201 (2013). 
[21] V. Fritsch et al, arXiv: 1301. 6062 

[22] E. D. Mun, S. L. Budko, C. Martin, H. Kim, M. A. 
Tanatar, J.-H. Park, T. Murphy, G. M. Schmiedeshoff, N. 
Dilley, R. Prozorov, and P. C. Canfield, arXiv:1211.063_6, 



